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Abstract. The general equation from previous work is specialized to a linear potential V{r) — —a + Fr 
acting in the space of spherically symmetric S wave functions. The fine and hyperfine interaction creates 
then a i-dependence in the effective potential energy equation and a position dependent mass m(r) in 
the effective kinetic energy of the associated Schrodinger equation. The results are compared with the 
available experimental and theoretical spectral data on the vr and p. Solving the eigenvalue problem within 
the analytically tractable Airy-function approach induces a certain amount of arbitrariness (fudge factors) . 
Despite of this, the agreement with experimental data is good and partially better than other calculations, 
including Godfrey and Isgur and Baldicchi and Prosperi |10| . The short comings of the present model 
can be removed easily in more elaborate work. 
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1 The S-state Hamiltonian 

For spherically symmetric S states the previously derived 
Hamiltonian reduces in Fourier approximation to [T]P| |^ 
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There are no more interactions than the central potential, 
the hyperfine, the kinetic, and the Darwin interaction, but 
also no less. For s-states the total spin squared is a good 
quantum number = [{ai + (T2)/2]^ = S{S + 1), thus 



0-1(72 = 2S{S +1) - 3 = 



-1, for 5 = 1, triplet, 
-3, for 5 = 0, singlet. 



(2) 



Because it is shorter, aia2 is kept explicit in the equations 
as an abbreviation for Eq.|21). Choosing a linear potential. 



V{r) = -a + Fi 



(3) 



with the force parameter F, often called string tension 
cr, the Hamiltonian ^ becomes a non-local Schrodinger 
equation with a ^-potential 
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Shaping notation, the Hamiltonian is written as 

H = ^ , , +Fr-a 
2mr[r) 

r \2 

The dimensionless 'coupling constant' is 

F 
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The spin-averaged potential energy for equal masses, 

r 
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has an attractive Coulomb potential. It has its origin in 
the Darwin term. The non locality of the Hamiltonian 
resides in the position dependent mass 
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To solve this Hamiltonian, one must go on a computer. 

The Hamiltonian in Eq.Q looks like a conventional 
instant form Hamiltonian as obtained by quantizing the 
system at equal usual time. But it must be emphasized 
that it continues to be a genuine front form or light cone 
Hamiltonian 0], derived from the latter by a series of exact 
unitary transformations ^E]. 



2 The model Hamiltonian and its parameters 

In this first round, I try to avoid to go on the computer as 
far as possible, by the following reason. The parameters in 
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Table 1. Model parameters in GeV. Note: [/*] = 1. 



fr fi 



ruh 



lOF 



2 1 0.4259 0.5553 1.8152 5.2505 1.1317 2.1454 



the theory must be determined from experiment, and this 
turns out as a non trivial, strongly non linear problem. In 
order to get a first and rough estimate, the Hamiltonian 
is simplified here until it has a form which is amenable to 
analytical solution. Therefore, all in-tractable terms in the 
above will be replaced here by mean values and related to 
the experimentally accessible mean square radius (r^) 'S^. 
In effect, the substitution 
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is the only true assumption in the present model. I con- 
sider thus the model Hamiltonian, 
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with the abbreviations 
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Its eigenvalues are 



En = —a + uj ^0 + uj r]n + ca\a2 ■, 



2mr 



(10) 

(11) 
(12) 
(13) 
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with the negative zeros of the Airy functions ^„ and r]n — 
^„ — ^0- A few ones are tabulated in App.^ The invariant 
mass squares 

Mn = (mi -t-m2)^ 

-I- 2(mi + 7712) ( - a 4- + ?7rtW + caia2) , (15) 

are then related to experiment. 

For equal masses mi = m2 = m, the model has the 3 
parameters m, F and a. One thus needs 3 empirical data 
to determine them. I choose: 
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The spectrum is labeled self explanatory by the flavor 
composition M„ = Mdu,tn or M„ = Mau,sn, for singlets or 
triplets, respectively. The triple chosen in Ea. 1)16(1 exposes 
a certain asymmetry. The excited p is chosen since its ex- 
perimental limit of error is very much smaller than the 
one for the corresponding tt state. Only its ground state 
mass is known very accurately, i.e. mT^+ — 139.57018 ± 



Table 2. Dependence on the fudge factors. 
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0.00035 MeV. In the present work only the first 4 digits 
are used. For equal masses, the above abbreviations be- 
come 
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The experiment defines 2 certainly positive differences: 
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A third one will be constructed by the observation that 



j^^-- i^^~-^^dji,sO — 4ma-4m^. Keeping in mind that 
^3 _ F"^ /fh, one can remove trivial kinematic factors and 
define 3 experimental quantities -B, C and D by 
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Substituting F — mC and ma ~ ' ™2 



m, gives 
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a quadratic equation with the solution 
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Having m, the F and a are then calculated from H19|l . 

With Airy functions, moments of different powers in 
r are somewhat difficult to evaluate. Therefore Gaussian 
weights are used, which give 



_8_ 
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In order to allow for corrections due to the true wave func- 
tions I introduce two fudge factors /* according to 



(22) 
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Table 3. S wave spectra in GeV for light unflavored mesons. 



Table 4. S wave spectra in GeV for strange mesons. 
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Hagiwara etal 7 , Zhou and Pauli 8_. 
^Godfrey and Isgur (9), *Baldicchi and Prosperi |1U| (a), 
"Could be a D state [TT|. 

Since all masons have about the same size by order of 
magnitude, these numbers are kept universal. The fudge 
factors are introduced here to account, in some global fash- 
ion, for the tremendous simplification introduced by re- 
placing Eq.lO with (|10|l . Some large scale variations of 
/* and /* are compiled in Table |21 The mass spectra in- 
cluding the ground states vary very little with the fudge 
factors. Any variations would show up the fastest for the 
high excitations. For this reason, the masses for n = 4 are 
included in the table. I do not understand this insensitiv- 
ity from a mathematical or numerical point of view. The 
major effect of /* is the ease by which one can change the 
quark mass. A value of /* ~ 40 leads to the 20 MeV for 
the quark mass quoted in JH]. In the present model, the 
values /* — 2 and /* = 1 are taken without seeking an 
optimum. 

In principle, one could determine the heavier quark 
masses analytically from the hyperfine splittings. The so 
obtained results are, however, not very reasonable, since 
experimental numbers are not sufficiently accurate. There- 
fore, I determine them numerically from Mus,sOj -^ue.so 
and Muc.so and compile them in Tabled The force pa- 
rameter F ^ 1100 MeV/fm is in line with currently used 
string tensions. 



3 Results and Discussion 

Unflavored light mesons. The results for the ir-p sys- 
tem are compiled in Table |21 The experimental points are 
taken from from Hagiwara et al [7]. It is no surprise that 
theory and experiment coincide for the tt"*", the p"*" and 
the p+(1450), because these data have been used to deter- 
mine the parameters. More surprizing is, that the model 
reproduces the huge mass of the excited pion within the 
error limit. This solves for me a long standing puzzle: 
Why is it that the ground state has a mass of 140 MeV, 



n Experiment^ Theory n Experiment^ Theory 



^So Singlets K+{us) 
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0.90^ 
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1.9598 
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^Hagiwara etal 0, "^Zhou and Pauli 

^Godfrey and Isgur 0, 

"To be confirmed; ^J^ not confirmed. 



while the first excitation with its 1300 MeV is different 
by an order of magnitude? — The answer is the usual 
one: Two large scales interfer destructively for the ground 
and constructively for the excited state. The two scales 
are the depth a ^ 1100 MeV and the string constant 
F - 1100 MeV/fm. 

The remaining three experimental masses of the n-p 
sector agree reasonably well with the calculation. There 
is no confirmed datum for the second excited p+ (3'^Si). 
The third excited p+ (4'^Si) deviates from experiment 
by 18 MeV, which is almost wihin the limits of error. 
The second excited (3^So) deviates from experiment 
by 107 MeV. The table includes also a comparison with 
other theoretical calculations. It includes the results from 
a recent oscillator model [S]. Their model is even simpler 
than the present one: it works with a hyperfine splitting, 
only, but suppresses the mechanism of a position depen- 
dent mass. Despite this, their results differ little from the 
present ones. I have included also the results from the pi- 
oneering work of Godfrey and Isgur |^ as a prototype of 
a phenomenological model, and from a recent advanced 
calculation by Baldicchi and Prosperi JO]. Note that ei- 
ther of these models have not much in common with the 
present one. The present model gives a good description 
for the pion, particularly the small mass of the physical 
pion is reproduced. All other potential models, including 
even Baldicchi and Prosperi, have the wellknown difficul- 
ties with that. I would have loved to compare also with 
Lattice Gauge Calculations but of course no data are avail- 
able for exited states, particularly not for such light sys- 
tems as the pion and the rho. Note that the effort on a 
pocket calculator is ridicuously small as compared to gi- 
gaflops years of calculations. 

Strange mesons. The S wave and K*^ spectra 
are given in Table0| The mass of the ground state of is 
used to determine the mass parameter nis- The excitations 
for the K (n^So) differ by only 60 and 3 MeV, respectively, 
and the spectrum for the K* (n^Si) by 27 and 64 MeV. 
Possibly, this could even be improved by playing with the 
fudge parameters, but in view of the experimental situ- 
ation, it is not done here. Except the ground states, the 
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Table 5. Ground state masses in GeV for heavy mesons. Table 6. The predicted S spectrum in GeV for heavy mesons. 



n Experiment^ Theory n Experiment^ Theory 



'So Singlet D"{uc) 



1.8645(5) 



1.8645 
1.9224^ 

1 oo3 



'So Singlet B+{ub) 



5.2790(5) 



5.2790 

5.2965^ 

5.31^ 



^So Singlet D~{sc) 



1.9685(6) 



1.9961 
2.0201^ 



'So Singlet B^'(s6) 



5.3696(24) 



5.3961 

5.3739^ 

5.35^ 



'So Singlet B+{cb) 



6.4(4) 



6.4914 



6.27^ 



'Si Triplet _D*"(mc) 



2.0067(5) 



1.9568 

2.0067^ 

2.04^ 



'Si Triplet B'+iub) 



5.3250(6) 



5.3082 

5.3250^ 

5.37^ 



'Si Triplet D^' {sc) 



2.1124(7) 



2.0665 

2.0655^ 

2.13^ 



'Si Triplet B;"(s6) 



5.4166(35) 



5.4185 

5.3885^ 

5.45^ 



■'Si Triplet Bt+{cb) 



6.4984 

6.3458^ 

6.34^ 



Hagiwara etal [3, ''Zhou and Pauh |S]. 
^Godfrey and Isgur 13, 
"To be confirmed; 'j^ not confirmed. 
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6.4733 
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6.4971 
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1.9961 
2.6878 
3.1427 
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2.0665 
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1 
2 
3 
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5.3961 
6.0349 
6.5114 
6.9052 


1 
2 
3 
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5.4185 
6.0549 
6.5299 
6.9227 


n^So 


Bt 


{cb) 


n^Si 


b:+ 
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1 
2 
3 
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6.4914 
6.9688 
7.3365 
7.6467 


1 
2 
3 
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6.4984 
6.9752 
7.3426 
7.6526 



experiments carry many ambiguities about the quantum 
number assignment for K and K* mesons. Both the first 
and the second excited state of K (2^So and 3^So) are 
not confirmed. Another unconfirmed resonance with mass 
1.629 ± 0.027 GeV lying between 2^80 and S^Sq was as- 
signed to be a singlet K . Apparently there is no position 
for it in the K spectrum if it is an S wave state. However, 
according to its mass and the present work, it might well 
be the first excited state of K* (2^So). 

Heavy mesons. The S wave uc, ub, sc, sb and cb 
meson spectra are given in Table |5l No excitations were 
observed for these mesons. — The mass of the ground 
state of is used to determine the mass parameter rric. 
No much data are available for D and D* mesons. The 
model prediction for I^Sq of D*'^ is smaller than experi- 
ment by 50 MeV. — The mass of the ground state of 
is used to determine the mass parameter mi,. The ground 
state of B*'^ agrees with the present model to within 17 
MeV. — No experimental values in the sc mesons are used 
to determine the model parameters. The model deviates 
from the available ground states by 27 and 45 MeV. — 
No data in sb mesons are used to determine the model 
parameters. The model deviates from the experiment by 
26 and 2 MeV, and is thus almost within the experimen- 
tal errors. — The mass of the ground state of (I^Sq) 
carries a large experimental error. Model and experiment 
agree. — The model prediction for the excited states are 
compiled in Table for easy reference. 



4 Conclusions 

The agreement between the present simple model and the 
experiment is excellent, with small but significant devi- 
ations. Perfect agreement has not been the goal of the 
present work. There must be room for a possible improve- 
ments by the 'true' equation jSJ- 

With the 4 mass parameters of the up/down, strange, 
charm and bottom quarks, the model has only 2 additional 
parameters for the linear potential. In principle, the fudge 
factors should be counted as parameter as well, but as 
seen above, the choice of the up/down mass and the fudge 
factors is strongly coupled. Thus, with 6 canonical param- 
eters the model exposes a reasonably good agreement with 
all 21 available data points. 

Note that renormalized gauge field theory has also 
4-f 1-1-1 parameters: The 4 fiavor quark masses, the strong 
coupling constant as, and the renormalization scale A. Of 
course, they can be mapped into each other PE]- 

Once one has determined the parameters in such a first 
guess, one should relax the model assumption, Eq.(|5J, and 
work with the full non local model, with a position depen- 
dent mass. For this one has to go back to the computer 
and perform the necessary fine tunings of the parameters. 
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A Solution to the linear potential 

Restricting to spherical symmetry (/ = 0), the Schrodinger 
equation for V{r) = Fr is 



p 

2mr 



1 1 

2mr r dr^ 



Fr 
Fr 



En 
En 



(23) 



<'(r) - 2mriFr - En)vn{r) = , 



In the last step w„(r) = r!f'„(r) was substituted. With the 
dimensionless variable 



(24) 



equation (|23|l is mapped into the differential equation for 
the Airy function mJ'(^) — ^Un(C) = 0- Eigenvalues are 
obtained by boundary conditions, in this case w„(0) = 0, 



En — 



2mr- 



= flUJ £.n 



(25) 



tf'„(r) = ^Ai{r[2mrF]i - Q . 
The ^n are the negative zeros of the Airy functions IH]: 



n 




Vn 







2.33811 








1 


4.08795 


1.74984 


1.74984 


2 


5.52056 


3.18245 


1.43261 


3 


6.78670 


4.44853 


1.26614 



(26) 



The table includes rjn ~ S,n ~ "^o- The table includes also 
the spacings Sn = £,n — Cn-ij which vary slowly with n. 



